Polar modes in multiferroic BaiMg^FeyiO^ hexagonal Y-type 
hexaferrite 



G.A. Komandin 1 , A.S. Prokhorov 1 ' 2 , V.I. Torgashev 3 , A. A. Bush 4 , E.S. Zhukova 1 ' 2 , and B.P. Gorshunov 1 ' 2 ^ 

1 A.M.Prokhorov Institute of General Physics, Russian Academy of Sciences, Vavilov str. 38, 119991 Moscow, Russia 

2 Moscow Institute of Physics and Technology (State University), Institutskii lane 9, 141700, Dolgoprudnyi, Russia 

3 Faculty of Physics, Southern Federal University, Zorge Str. 5, 344090 Rostov-on-Don, Russia 

4 Moscow Institute of Radiotechnics, Electronics and Automation, 117464, Moscow, Russia 

the date of receipt and acceptance should be inserted later 

Abstract. Spectra of transmission and reflection coefficients of single crystalline F-type hexaferrite 
Ba2Mg^F e\iOi2 are measured at Terahertz and infrared frequencies, 3 — 4500 cm" 1 , and at temperatures 
10 to 300 K, for polarization E _!_ c. Oscillator parameters of nineteen phonon lines allowed by the R3m 
crystal lattice symmetry are determined. Absorption lines assignment to vibrational modes of the lattice 
structural fragments is performed. Below the temperatures of 195 K and 50 K, corresponding to the zero 
field phase transitions to the proper screw and to the longitudinal-conical spin states, respectively, new 
absorption lines are discovered whose origin is assigned to the electric-dipole active magnetic excitations. 
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1 Introduction 

After the discovery of a huge magneto-electric effect in TbMnO^ [1 , much attention is being paid to a class of 
ferroelectrics with cycloidal spin order |2I3| . In such materials the application of relatively small magnetic field (few 
Tesla) can induce various phases including those revealing ferroelectricity. Model objects for studying magnetic and 
ferroelectric effects are hexaferrites, typical representatives are compounds of the M-type, like AFe^Oig, Y-type, like 
A 2 B 2 Fei2022, W-type, like AMe2Feie027, Z-type, like A-iMe^Fe^O^x, X-type, like A 2 -Me2i' 1 e28046, an d U-type, 
a E-mail: gorshunov@ran.gpi.ru 
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like A^Me2Fes§0§Q\ here A = Sr, Pb, Ba; M = Zn, Fe, Co, Mg, Mn 0]. Of special interest are Y-type materials 
in which unique multiferroic properties have recently been discovered [516171819110111112113114115116] . In particular, in 
B0L2M giF enOii (BMFO), magnetoelectric switching can be realized by application of rather small magnetic field 
B (the value of B can be as low as ±30 mT [6 ), while Bao^Sri^Zn2(Fei- x Al x )i2022 can be efficiently tailored by 
the AZ-substitution (at x — 0.08 the critical magnetic field for switching the electric polarization is reduced down to 
~ 1 mT |12j). The phase diagrams of these compounds contain several magnetically ordered phases and some of them 
are also ferroelectric. The crystal structure of Y-type hexaferrites is rather complicated. For instance, the unit cell of 
B0L2M2F ei2C>22 (M = Zn, Fe, Co, Mg, Mn) contains 18 atomic layers with (cchhhh)^ stacking sequence along the 
c-direction |17ll8j . as indicated in figure 1. 

Magnetic structure of Y- hexaferrites is under intensive study for already more than 50 years [13115116119120121122123124"] . 
In |19[ . an existence of a collinear-spin ferrimagnetic order was suggested. Later, an idea about helical spin ordering in 
(Sr, Ba)2Zn2F e\2 022 was put forward [5D] and experimentally confirmed in [151161211 22 23 24 J. According to [TS], the 
magnetic structure in BC12M g2F 'evi022 is proper screw below 195 K at B — (figure lb) and collinear at higher tem- 
peratures (figure la). Application of magnetic field can generate a series of magnetic phase transitions [617110111113] . In 
Ba,2Mg2Fei2022 there exists a longitudinal conical spin ordered state below 50 K (figure lc) 13 . Here, ferroelectric 
phases with a quasitransverse spin ordering [6 13 (figure Id) can be induced by applying magnetic field of magnitude 
much smaller than in BaQ^Sr\^Zn2Fe\2022 |5I8) . It turned out that the phase boundaries can be effectively tuned 
by heat treating [5] as well as by ^/-substitution into octahedral Fe sublattices 12 . It should be noted that recent 
neutron diffraction studies |10I13[ revealed much more complicated magnetic phase diagram of Ba2Mg2Feyi022 than 
was reported previously [6] . 

Since hexaferrites with magnetically ordered long periodic structures reveal strong magnetoelectric effects at ele- 
vated temperatures (above 77 K) and in relatively small magnetic fields, they are promising for application in magne- 
toelectric systems. Presently, the mechanisms of competition between magnetic and electric interactions in hexaferrites 
remain unclear. That is why the aim of this work was to further study the mechanisms of formation of magnetic and fer- 
roelectric phases in Ba2Mg2Feyi022 by using the Terahertz-infrared (THz-IR) spectroscopy. Indeed, finite interaction 
between the two subsystems, spins and phonons, would necessarily lead to a noticeable transformation in the phonon 
subsystem (IR range) in the materials and these changes can be used for getting more information about magnetoelec- 
tric properties of the compounds. Besides, recently, Kida et al. have reported on the observation in BaiMgiFey^Oii 
of electric-dipole active magnetic excitations [TT] . called electromagnons. The excitations were detected at THz fre- 
quencies in longitudinal-conical spin state and revealed huge changes in the intensity under application of magnetic 



G.A. Komandin et al.: Polar modes in multiferroic BaiMg-zFevzOyi hexagonal Y-type hexaferrite 



CD 01 O Fe © Mg(Fe) 

© 03 j Fe(a) © Ba 

® 04 O Fe(b) 

© 05 ® Fe(c) 

© 08 • Fe(Mg) , j 



(b) 

k„//[001] 




(c) 



<2> 



a B=0 



B=0 



(d) 

k 



Fig. 1. (a) - Unit (hexagonal) and primitive (rhombohedral) cells of Ba-iMgiFenOii hexaferrite. Coordinating 
polyhedrons of some cations in different cell sites are shown. T and S stand for two kinds of magnetic sublattices with 
large and small magnetic moments, (b) - (d) Schematic illustration of helix spin ordering in different magnetically 
ordered phases, (b), (c) and (d) correspond to the proper screw (50 < T < 195 K), longitudinal conical (T < 50 A") 
and quasitransverse (T < 195 K, B ^ 0) spin orderings, respectively. Spins orientation in the ferrimagnetic phase is 
indicated by arrows. 



field less than 0.2 T (along [001]). On the other hand, Kamba et al. report [Tl] that in Bao.^Sri ^ZniFenOii "no 
evidence for ferroelectricity was observed either on cooling (at B — 0) or after application of external magnetic field up 
to 0.3 T at 300 K". To understand the reason of these discrepancies was also one of our aims. 
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2 Experimental details 

Single crystalline samples of 3a,iMg<iFzviOi<2, were grown by a flux method as described in |16j . For the infrared 
and THz measurements plane-parallel samples were prepared of thickness about 0.5 mm and size about 5x5 mm 2 , 
with the c axis perpendicular to the surface of the samples. In this geometry, only the E u modes can be registered 
by optical measurements, corresponding to the situation when the .E-vector of the probing radiation is perpendicular 
to the c-axis. The measurements were performed on two spectrometers. For the low frequency range, 15 — 35 cm -1 , 
a quasioptical THz spectrometer based on backward-wave oscillators was used. With this instrument by measuring 
complex transmission coefficient spectra Tr{u), the values of real e' and imaginary e" parts of dielectric permittivity 
are directly obtained [25] . A Fourier-transform spectrometer Bruker IFS 113V was used to cover higher frequency range 
up to 4500 cm^ 1 . At these frequencies the reflection coefficient spectra R{v) were measured. The spectra obtained on 
the two spectrometers, Tr(y), R(is), e'(^) and e"(v), were merged and analyzed by a least square fitting procedure. 
The absorption lines were modeled with Lorentzians and coupled Lorentzians with the complex dielectric permittivity 
given by the following expressions. For Lorentzians 

n A 2 

Ui A — V 1 + IVY; 
3 = 1 J J 

where £oo is the high frequency dielectric constant, Aej is the dielectric contribution of the absorption line, Vj is its 
eigenfrequency and the damping constant. For the coupled Lorentzians the corresponding expression is written as 

^ ^ _ S 2 (v 2 2 - V 2 + W7 2 ) + 52 (vi 2 -V 2 + ivil) - 2 y/5152 (g + wS) 

{vi 2 — v 2 + iv~fi) (v2 2 — v 2 + W72) — (a + iv8) 2 
where j = 1, 2, Sj = AejV 2 is the oscillator strength of j-th Lorentzian with the eigenfrequency v$\ a is the real and 
d the imaginary coupling constants. 



3 Results and discussion 

Examples of the spectra of the -Ba2M(? 2 Fei 2 022 sample obtained at 10 K are presented in figure 2. Dots show the 
spectra got on the THz and Fourier spectrometers and the lines represent the results of the least-square processing 
of the spectra based on the expressions [T] and [5] At the lowest frequencies, an oscillating behavior in the spectra is 
observed that comes from interference of the radiation within the plane-parallel samples. Processing the oscillations 
allows to determine directly the dielectric parameters of the material, the real and the imaginary parts of the dielectric 
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constant, of the refractive index, etc. [25 J. Figure 3 shows examples of the e' and e" spectra of -Bo^Mf^Fe^C^- Dots 
correspond to the values of e' and e" directly obtained on the THz spectrometer [25] while lines show the results of 
the fitting of the infrared spectra. Parameters of the absorption lines are summarized in [T] 

3.1 Spectral analysis 

We have performed a factor-group analysis of possible phonon modes in Ba2Mg2Fe\2 022, with the results summarized 
in table 2. According to this table, 34 modes can be observed in the IR-spectra corresponding to the irreducible 
representations 15^2^ \IR ■ z] + 19E U [IR : x,y] of the space group Rim. Nondegenerate ^^-modes can be excited 
by electric field E \\ c and doubly degenerate F„-modes - by the .E-field polarized perpendicular to the c-axis. 
To identify the excitations observed in our experiment, we compare spectra of Bc^M ^-Fe^C^ with the ones of 
magnetoplumbite and spinels [27 28 29 30 31 32 33]. We consider two types of tetrahedra (M6 and M9 in table 2). An 
ideal FeO^ tetrahedron with the Td symmetry has four fundamental vibration modes |34j: v\(A{), ^(F), 1/3 (Fa), and 
^4(^2), of which only V3 and V4 are IR-active. Due to the low local symmetry of Be^-M^^e^C^, also the V2 modes 
are active in optical spectra. We assume that the lines seen at 539 and 580 cm -1 are caused by the stretching V3 
vibrations and that the deformational modes of the tetrahedra [FejMg)0^ and (Mg/Fe)Oi are located within a 
complicated absorption structure in the range around 419 cm" 1 . The modes V2 have small dipole moment and should 
be observable in the range 300 — 350 cm -1 . The v\ mode is not IR-active. The complicated shape of the V3 and V4 
bands can be caused by a LO-TO splitting of the polar modes or by a certain disorder of cations (Fe/Mg and Mg/Fe) 
at the (6c) sites within the tetrahedron leading to different eigenfrequencies. 

An undistorted FeOs octahedron has a Oh symmetry and creates 6 fundamental modes [34j : Vi(Ai g ), v^iEg), 
1/3 (.Flu), v i{Fiv)i l/ 5{F2 g ), and v^iFiu)- Vibrations v\, vi, and v§ should be Raman active, V3 and - IR-active and 
is expected to be silent. Since the FeOs octahedra form chains and thus are not independent, unambiguous assignment 
of the observed absorption bands to octahedral v\ — Vf, modes is problematic. That is why only tentative assignment of 
these bands is given in table 1. Still, it is reasonable to assume that the bands at 240 — 300 cm -1 , 350 — 400 cm" 1 and 
442 m _1 correspond to the 1^3 and 1^4 vibrations of the octahedra. The bands at 150 — 210 cm~ l are mainly coming 
from the bending vibrations of the concatenated O — Fe — O — Fe — O bonds of adjoint octahedra and could also 
be assigned to the V3 and ^4 molecular vibrations of the octahedron. In 30J it is assumed that the absorption bands 
located close to 150 cm -1 and 175 cm" 1 in the spectra of M-type hexaferrite BaFe^Oig are caused by vibrations 
of the bipyramid. The structure of Ba,2Mg2Fei2 022, however, does not contain bipyramids and thus our assignment 
of the bands in this range to the bending vibrations looks more correct. Finally, the line around 92 cm -1 should be 
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Fig. 2. Spectra (dots) of reflection (a) and complex transmission (amplitude and phase shift) (b) [25] of the 
BaiMg2Fe\iOii sample measured at T = 10 K for polarization E _L c, and their model processing (lines, see 
text). Filled symbols in panel (a) are obtained by processing the THz transmission coefficient and phase shift spectra 
of the type shown in panel (b) . 



assigned to the Ba — O bonds vibration. We note that in BaFe^Oig this vibration has slightly larger frequency, of 
100 - 102 cm' 1 and 123 - 125 cm' 1 , according to [28130131] . 



3.2 Low temperature efFects 

We have performed measurements of the THz-IR spectra of Ba>2M ^^viOii down to liquid helium temperatures. 
Figure 3 shows corresponding spectra of e' and e" at room temperature and at 10 K . Two features have to be pointed 
out. 
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Table 1. Parameters of modes detected in the spectra of single crystalline Ba^Mg^FenOii for polarization E _L c 
at T = 10 K: frequency position i/j, damping jj, dielectric contribution Aej; a and S characterize the modes coupling 
between the modes. The last column shows assignments of the modes to the type of displacement of the structure- 
forming atoms. The high-frequency dielectric permittivity is = 4.52. 



Oscillator type 


Vj , cm 1 


jj , cm 1 




a 


S 


Assignment, ions involved 


coupled Lorentzians 1 


29.08 


20.52 


0.36 






magnon-phonon continuum 


56.78 


41.32 


0.36 




-0.86 


Lorentzian 2 


93.21 


1.62 


0.67 







Ba-O 


Lorentzian 3 


130.41 


8.27 


0.20 






Fe-O-Fe, {FeOe) 


Lorentzian 4 


140.65 


7.41 


0.23 






optic magnon 


Lorentzian 5 


152.93 


9.11 


0.30 






Fe-O-Fe, {FeOe) 


Lorentzian 6 


167.45 


9.11 


0.21 






Fe-O-Fe, {FeOe) 


Lorentzian 7 


176.29 


7.96 


0.75 






Fe-O - Fe, (FeOa) 


Lorentzian 8 


211.54 


17.54 


1.40 






Fe-O-Fe, (FeO G ) 


Lorentzian 9 


233.52 


19.55 


1.21 






optic magnon 


coupled Lorentzians 10 


243.14 


8.93 


0.45 







Fe-O, {FeOe), bend 


288.47 


11.10 


1.13 




-43.38 


Fe - O, {FeOe), bend 


coupled Lorentzians 11 


301.30 


32.57 


2.50 







Fe-O, {FeOe), bend 


325.10 


6.54 


0.87 




6.73 


U2 {FeO^), bend 


coupled Lorentzians 12 


417.59 


13.21 


0.75 







v A {Fe/Mg)0 4 , bend 


431.33 


11.58 


0.12 




-13.11 


v A {Mg/Fe)0 4 , bend 


coupled Lorentzians 13 


353.49 


28.33 


0.09 







Fe - O, {FeOe), bend 


403.86 


29.65 


0.80 




38.23 


Fe-O, {FeOe), bend 


coupled Lorentzians 14 


447.81 


30.79 


0.29 







Fe - O, {FeOe), stretch 


520.16 


24.38 


0.16 




-0.26 


Fe - O, {FeOe), stretch 


coupled Lorentzians 15 


538.47 


22.55 


0.56 







i>3 {Fe/Mg)0 4 , stretch 


582.43 


20.26 


0.39 




-6.00 


v 3 {Mg/Fe)0 4 , stretch 



22 



^e(0) = e ao + Y, Ae J = 18 - 32 
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Fig. 3. Spectra of real (a) and imaginary (b) parts of dielectric permittivity of Ba2-M c/2-F zviOii single crystal at 
300 K and 10 K, for polarization E _L c. Dots at low frequencies (below 40 cm -1 ) show the data obtained from 
transmission coefficient spectra of a plane-parallel sample, see text. 

Firstly, the low frequency (around 15 cm -1 ) value of e' does not change significantly while cooling down, with 
the total dielectric contribution of all higher frequency absorption bands of At w 18. This means that there are no 
significant structural changes in the crystal lattice of Ba2Mg2Fe\?,Oii while cooling down except trivial narrowing of 
the bands and a certain shift of the line around 130 cmT 1 (figure 4b). 

Secondly, we were able to register for the first time several new absorption lines appearing below 195 K and below 
50 K. These two specific temperatures correspond to the change in the magnetic structure of BaiMg^FenOii |13j . 
The new lines are indicated by arrows in figure 3 and shown in more details on expanded scales in figure 4. 

At the temperature of 195 K, the configuration of the internal magnetic structure of Ba2M 32 F ei2022 changes from 
collinear ferrimagnetic at T > 195 K (figure la) to an incommensurate proper screw spin-ordered at 50 < T < 195 K 
[13] with the wavevector ko || [001]. This magnetic transition leads to an appearance of three new lines at 85 cm -1 , 
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Table 2. Factor-group analysis of mechanical representation for Ba^M g2f ^e-vzOyx hexaferrite. In the first column, the 
number at the atomic symbol represents the sublattice type. 

Atom Wyck.pos. Site symmetry Vibration representation Atomic environment 



01 


18h 


m[C s (6)] 


2A lg + Am + A 2g 


+ 2A 2u 


+ 3E g - 


^3E U 


coplanar triangle Fez 


Fe2 


18h 


m[C s (6)] 


2A lg + Am + A 2g 


+ 2A 2u 


+ 3E g - 


h 3E U 


octahedron C?6 


03 


18h 


m[C s {6)] 


2A lg + Am + A 2g 


+ 2A 2u 


+ 3E g - 


h3E u 


non-coplanar triangle Fe 3 


04 


18h 


m[C s {6)] 


2A lg + Am + A 2g 


+ 2A 2u 


+ 3E g - 


h3E u 


tetrahedron MgFe 3 


05 


6c 


3m[C 3 „(2)] 


Alg + A 2u +Eg + 


E u 






tetrahedron Fa 


M6 


6c 


3m[C 3v (2)] 


Alg + A 2u +E g + 


E u 






tetrahedron O4 


Ba7 


6c 


3m[C 3 „(2)] 


Alg + A 2u +Eg + 


E u 






9-vertex polyhedron O9 


08 


6c 


3m[C 3 „(2)] 


Alg + A 2u +Eg + 


E u 






tetrahedron MgFes 


M9 


6c 


3m[C 3v (2)] 


Alg + A 2u +Eg + 


E u 






tetrahedron O4 


FelO 


6c 


3m[C 3v (2)] 


Alg + A 2u +E g + 


Eu 






octahedron 0§ 


Fell 


3b 


-3m[D 3d (l)] 


A 2u + E u 








octahedron Oa 


Fel2 


3a 


-3m[£> M (l)] 


A 2u + E u 








octahedron 0% 



M6 = 0.75Fe + 0.25Mp; M9 = 0.75Mg + 0.25Fe 



Fmech = 14Ai s [Raman : (a xx + a yy , azz)] + 4A lu (silent) + AA 2g [R z ] + l6A 2u (IR : z)+ 

+ 18E g [(Rz, Rz); Raman : (a xx - a yy , a xy ); (a xz , ayz)] + 20 E U [IR : x,y]. r acous u c = A 2u + E u . 

r me ch = UA lg [Raman : (a xx + a yy ,azz)} + 4 Am (silent) + 4:A 2g [R z ] + 15A 2u (IR : z) + 

+ 18E g [(Rz, Rz); Raman : (a xx — a yy , a xy ); (a xz ,ayz)] + 19E U [IR : x,y). 



140 cm -1 and 232 cmT 1 (figures 3, 4). We believe that these modes have magnetic origin and are related to an 
excitation of optically active magnetic polaritons 35 55]. To check our suggestion, optical measurements in magnetic 
field are in progress. 

Another excitation was discovered in B-C12M e^C^ at the lowest frequencies, below 40 — 50 cm~ 1 . As shown in 
figure 4c, there is a broad wing below 80 cmT 1 in the e" spectrum at 100 K that transforms into a more pronounced 
absorption band at 10 K. In [11] by using the technique of THz time-domain spectroscopy, an absorption band has 
been detected in Ba,2M enO-zi at 22 cm -1 below 50 K (where the magnetic structure has changed to a longitudinal 
conical [13]) and has been assigned to an electromagnon excitation. Although this line was observed in the geometry 
E || c, we assume that the transition around 50 K, that leads to an activation of an electromagnon for E || c, 
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Fig. 4. Spectra of imaginary part of dielectric permittivity of BaiMg2Fe\iOii single crystal at different temperatures, 
for polarization E _L c. The spectra demonstrate appearance of new absorption lines below the magnetic phase 
transition temperature T rs 195 K. 

is accompanied by certain distortions in the crystal lattice of Ba2MgiFe\202i which break the selection rules and 
enable excitation of the mode also for polarization E _L c. 

According to the spin currents model |37j (or the inverse Dzyaloshinskii-Moriya model), the non-collinear spin 
ordering (figure Id) can produce a local polarization p^j = Aetj x (Si x Sj), that is vanishing on a macroscopic scale. 
Here the constant A depends mainly on the spin-orbit and the exchange interactions. The neighboring spins Si and 
Sj, that connect sites i and j, are aligned along the unit vector eij = feg/ | feo |- Within this model, in order to create 
the ferroelectric polarization, the alignment of the cone axis should deviate from the direction feo- (It is important 
that such changes should happen, on the microstructural level, to produce the value of A not equal to zero). The 
Ba-zM ziiO-zi crystal is composed of an alternating stacking of the T and S blocks along the c axis, as depicted 
in figure la. Within these blocks, the magnetic moments on Fe ions are collinearly aligned and ferrimagnetically 
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Fig. 5. Displacements corresponding to the bending mode and superexchange interactions in BaiMgiFe\iOi2- 

coupled. Along with that, the magnetic moments between the adjacent T and S blocks align noncollinearly (figure 
5). The resultant magnetic structure of BaiMg^Fey^O^i in the absence of magnetic field is a noncollinear screw 
spiral magnetic structure. This can lead to local distortions of magnetic order and to nonzero values of A. When the 
temperature is lowered, the exchange interactions change their amplitudes (in particular, on the boundaries between 
the T and the S blocks), probably due to spin-phonon interaction with the deformation Fe — O — Fe modes involved 
(figure 5a) . It is this fact that explains the evolution of shape of the spectral response presented in figure 4. 



4 Summary 

By using the terahertz frequency-domain and the infrared Fourier spectrometers, the spectra of the real and the 
imaginary parts of the dielectric permittivity have been measured of a single crystalline Y-type hexaferrite 
Ba2Mg 2 Fei2 022, at frequencies 3 cm^ 1 to 4000 cm -1 in the temperature interval from 5 K to 300 K and for 
polarization E _L c. 

New absorption lines are discovered below the temperatures of 195 K and 50 K, where the internal magnetic field 
configuration in Ba,2Mg2Fevi022 changes. The origin of these lines is assigned to optical polaritons. Their spectral 
weight is nonzero in the phase with longitudinal spin polarization (T < 50 K) testifying the influence of the spin-lattice 
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interaction (that involves the deformational Fe — O — Fe modes) on the superexchange interaction at the junction 
between the S and the T structural blocks. 

The parameters of all nineteen (allowed by the i?3m symmetry) phonon E u modes are determined. Their assignment 
to the vibrational modes of the crystal lattice fragments is performed. 

It is found that the total contribution to the static dielectric permittivity of the polar phonon modes is practically 
independent on the temperature in the range 10 K — 300 K, implying no significant structural changes in the crystal 
lattice of BaiMgiFe\iOii- 
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